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氢化 物 对 钳 拉 伸 性 能 影响 的 分 子 动力 学 人 研究 


XE HEr KE 
(华北 电力 大 学 核 科 学 与 工程 学 院 北京 102206) 


摘要 氧化 物 是 错 合 金 包 壳 管 在 核电 厂 正常 运行 过 程 中 与 一 回路 冷却 剂 发 生 错 水 反应 而 产生 的 常见 缺陷 。 本 文 利用 分 
子 动力 学 方法 , 采用 COMB3 势 函 数 , 构建 含 氢 化 物 的 钻 基 模型 进行 单 轴 拉 伸 模拟 , 探究 了 気 化物 密度 対生 力学 性 能 
的 影响 。 研 究 结果 表明 ， 当 氢化 物 密度 在 0~1078 ugg 时 ， 随 着 氧化 物 密度 的 增加 ， 届 服 强度 、 应 变 和 杨 氏 模 量 降低 。 
在 弹性 阶段 ， 氨 化物 密 度 的 增加 使 应 力 集中 区 域 增 大 ， 有 利于 位 错 形 核 ， 在 塑性 变形 阶段 ， 随 着 氢化 物 密度 的 增 大 ， 
初始 位 错 更 倾向 于 在 氢化 物 周 围 扩 展 。 当 氢化 物 密度 在 1078 ~ 2311 ug/g 時 , 随 気化 物 密 度 的 増加 , 屈服 強度 、 成 変 
和 杨 氏 模 量 升 高 ， 这 是 由 于 氧化 物 密度 较 高 时 产生 了 大 量 位 错 并 造成 位 错 塞 积 。 

关键 词 ” 错 合金 ， 氧 化物 ， 单 轴 拉 伸 ， 分 子 动力 守 
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Abstract [Background]: Hydride is a common defect caused by the reaction of zirconium water with primary coolant in 
the normal operation of a nuclear power plant. [Purpose]: In order to investigate the effect of hydride density on mechanical 
properties of zirconium. [Methods]: Through the molecular dynamics simulation method and COMB3 potential function, we 
constructed the zirconium base model containing hydride and uniaxial stretching was carried out. [Results]: The results show 
that the yield strength, strain and Young's modulus decrease with the increase of hydride density when the hydride density 
ranges from 0 to 1078ug/g. When the hydride density is between 1078 and 2311pg/g, the yield strength, strain and Young's 
modulus increase with the increase of hydride density. [Conclusions]: When the hydride density ranges from 0 to 1078yug/g, 
in the elastic stage, the increase of hydride density increases the stress concentration area, which is conducive to dislocation 
nucleation. In the plastic deformation stage, with the increase of hydride density, the initial dislocation is more inclined to 
expand around the hydride. When the hydride density is in the range of 1078 ~ 2311yg/g, a large number of dislocations are 
generated due to the high hydride density, resulting in dislocation plugging. 
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钳 合 金具 有 优异 的 高 温 力学 性 能 、 良 好 的 耐 腐蚀 性 能 和 低热 中 子 吸收 截面 等 特点 巾 ， 作 为 核燃料 包 壳 
管 的 重要 结构 材料 ， 已 被 广泛 用 于 轻 水 反应 堆 核 电厂 中 Pg。 在 使 用 期 间 的 一 个 重要 问题 是 ， HARES 
回路 冷却 剂 会 发 生 钳 水 反应 产生 氧化 物 沉淀 09I9， 氧 化 物 沉淀 相 的 产生 通常 会 影响 钳 合 金 的 力学 性 能 10- 急 。 
近 些 年 ， 许 多 学 者 利用 分 子 动力 学 (Molecular simulation; MD) 方 法 模 撤 研 究 気化 物 対 筐 力学 性 能 当 
影响 。A Lloyd 等 人 03] 采 用 MEAM 势 函数 研究 了 氧 浓度 对 Zr 力学 性 能 的 影响 ， 结 果 表 明 : H 的 浓度 越 
高 ， 材 料 永久 变形 所 需 的 极限 抗 拉 强度 和 断裂 应 变 就 越 低 。Ghaffarian 等 人 04 采 用 MEAM 势 函数 研究 了 拉 
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采用 MEAM Hf 
的 增加 有 助 于 位 错 增殖 。 最 近 ，Mark J Noordhoek 55 A65 H 
并 验证 可 正确 计算 Zr 体系 力学 性 能 。 

值 符合 较 好 。 
JT HU ES IK SINE PH 


由 于 氧 原子 
势 函 数 模 拟 


在 本 


预测 了 6 - ZrH1.0-2.0 的 形成 能 、 


m Zr 中 嵌入 氧化 物 的 变形 机 制 ， 阐 明了 氧化 物 析出 对 锋 合 金 延展 性 能 的 影响 。Poulami 等 人 [5] 


层 错 能 与 DFT 计算 相 匹配 ， 


究 中 ， 继 Zhang 等 人 0 验证 COMB3 势能 够 


COMB3 Z5: IR 


1 ”计算 方法 


运用 LAMMPS 分 子 动力 学 模拟 程序 构建 模型 ， 模 型 大 小 为 10.3 nmx10.1 nmx20.6 nm 的 立方 体 ， 整 


个 模型 内 包含 92160 个 Zr 原子 的 HCP 单 晶 ， 如 图 
唱和 向 ， 并 且 在 三 个 方向 上 均 采 用 周期 性 边 
的 COMB3 AKA EAMHAIN, KA 


[0001] jc p 


Noordhoek 5& A 161 H 


唱 格 参 数 和 体 模 量 , 5 DFT ii 
恰当 


究 氧 化 物 密度 对 独力 学 性 能 的 影响 ， 从 微观 原子 层面 古 


E 氧 化 物 后 ， 进 一 步 利 用 


究 気 化物 対 筑 力 学 性 能 


1 所 示 。X、Y、Z 


界 条 件 ， 以 减 小 边界 效应 的 影响 ， 势 函数 采用 


究 了 氧 浓度 对 错 的 作用 机 理 ， 发 现 随 着 H 含量 的 增加 ， 材 料 断 裂 时 的 应 变 值 减 小 ， 这 是 
Hf Zr 的 COMB3 经 验 势 ， 运 用 此 
Zhang 等 人 0 用 COMB3 势 


Ei 
影响。 


分 別 沿 [1210]、[1010]、 


HS- ZrH1.s 相 的 立方 结构 ,构建 了 半径 为 15 


A， 高 度 为 4.785 A 的 氧化 物 圆 盘 ， 唱 格 参数 为 ab, c = 4785 A (Ha, B, y= 90.0?) ， 晶 体 类 型 为 上 


CC 结构 ， 氨 化 物 在 Zr 基体 
间距 为 15 A， 每 添加 一 


初始 晶 胞 为 标准 晶 格 参数 及 结构 9H， 且 需要 在 Zr 体系 


等 体积 的 6- 


物 和 5 层 氧化 物 的 示意 图 。 


一 


Ex 


ZrHis 沉淀 。 将 6- ZrHis 沉淀 与 Zr 


计算 第 


即 总 时 间 50 ps， 


应 变速 率 沿 


1000 步 输出 原子 相应 的 坐标 信息 、 
数据 利用 Origin 绘制 成 应 力 应 变 
设 应 变 为 2 ， 应 力 为 0 ，E 为 杨 氏 模 量 ， 


伸 方 向 的 杨 


步 采 用 NPT 系 综 , 使用 COMB3 势 函数 在 
并 使 用 
[0001] 方 向 进行 单 轴 拉 人 


ETX Y, ZÝ 


按 (0001)w/{111}s 的 取向 关系 放置 。 不 同 


戻 数 (1-15 层 ) 的 気化 物 同 距 均 勾 且 


戻 気 化物 , 模型 気化 物 提 高 154 ug/g。 表 1 为 氧化 物 密度 与 放置 氧化 物 层 数 的 对 


300K 的 真空 环境 


心切 出 一 个 盘 形 体积 ， 
切割 的 空心 体积 相 匹配 。 图 2 为 在 Zr 基体 


并 在 其 位 置 插入 一 个 
放置 1 层 氧化 


31 500000 步 ( 步 长 为 0.1 fs), 
ERE (CG) 方法 使 系统 达到 能 量 最 小 化 。 第 二 步 在 温度 为 300 K, 以 1010s 
， 总 步 长 为 100000 步 ( 步 长 为 0.5 fs) ， 即 50 ps， 为 50% 应 变 。 期 间 


向 上 的 力 、 应 变 及 应 力 、 原 子 势能 ， 将 运算 所 得 的 


| 线 图 ， 并 通过 对 应 力 


应 变 


氏 模 量 。 


o= Ee 
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KHE 
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小 二 乘法 进行 线性 拟 合 ， 得 到 杨 氏 模 量 计算 


1 


式 (2) 


则 方程 为 


公式 为 


的 E 和 分别 为 应 变 和 应 力 样本 的 平均 1 


，D 为 样本 个 数 。 


采用 


适应 


构 变 化 并 输出 位 错 线 长 度 统计 。 借 助 可 视 化 软 从 


(Color coding) 进行 应 力 云图 可 视 化 分 析 。 


近邻 分 析 方 法 并 输出 结构 数量 百分比 、 位 错 分 析 DXA 来 观察 拉 人 
FOVITO 显示 模拟 过 程 的 微观 结构 3 
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1 线 中 线性 阶段 的 斜率 进行 拟 合 ， 来 得 到 该 拉 


(1) 


(2) 


变形 过 程 中 的 微观 结 


通过 原子 筛选 着 色 


刘 潇 雅 等 : 気化 物 対 笛 拉 伸 性 能 影 画 的 分 子 効 力学 研究 
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图 1 HCP 结构 Zr 基体 (a) 晶 向 及 XYZ 轴 方向 ，(b) 初始 模型 
Fig.l HCP structure Zr matrix (a) crystal direction and XYZ axis direction, (b) initial model 


表 1 不 同 氢化 物 密度 对 应 氢化 物 层 数 


Table 1 Different hydride densities correspond to the number of hydride layers 


氧化 物 层 数 

Number of hydride 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
layers 

気化 物 密度 

Hydride density/ 154 308 462 616 770 924 1078 1232 1386 1540 1695 1849 2003 2157 2311 
ugg? 


图 2 不 同 密度 的 氧化 物 放 置 在 Zr 基体 中 (a) 一 层 5 - 気 化物 , (b) 五 层 5 - 気 化物 
Fig.2 Different densities of hydride placed in Zr matrix (a) a layer of hydride, (b) five layers of hydride 
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2 ”结果 分 析 与 讨论 
2.1 ”应 力 应 变 曲线 及 力学 参数 


3 为 不 同 氧化 物 密度 下 的 Zr 基体 所 对 应 的 应 力 -应 变 曲线 。 通 过 图 3 可 以 发 现 所 有 模型 的 应 力 -应 变 
曲线 具有 相似 的 趋势 且 应 力 随 着 应 变 的 增加 而 增加 ,直到 达到 一 个 峰值 应 力 , 而 峰值 应 力 与 塑性 变形 有 关 ， 
并 且 把 应 力 -应 变 曲线 分 为 弹性 阶段 和 塑性 阶段 两 部 分 。 到 达 峰 值 应 力 以 后 ， 应 力 突然 下 降 到 一 个 平稳 值 ， 
随 着 应 变 的 继续 增加 ， 应 力 围绕 这 个 平稳 值 上 下 波动 。 

图 4 为 不 同 氧化 物 密度 下 的 屈服 应 力 和 杨 氏 模 量变 化 曲线 ， 在 氢化 物 密度 为 1078 ugg 时， 曲线 上 出 
现 拐 点 。 在 氧化 物 密度 于 0-1078 ug/g 时 ， 屈 服 应 力 、 屈 服 应 变 和 杨 氏 模 量 都 随 着 氧化 物 密度 的 增加 而 降 
低 。 氧 化 物 密度 为 1078 pg/g 时 , 模型 的 屈服 应 力 下 降 至 最 低 值 7.69 GPa, 比 纯 Zr 模型 届 服 应 力 下 降 42.2296; 
屈服 应 变 下 降 至 最 低 值 0.0895， 比 纯 Zr 模型 下降 39.34%; 杨 氏 模 量 下 降 至 最 低 值 112.18 GPa， 比 纯 Zr 模 
型 下降 8.94 %。 氢 化 物 密度 为 1078-2311 ug/g 时 ， 模 型 的 届 服 应 力 、 届 服 应 变 和 杨 氏 模 量 随 着 氢化 物 密度 
的 增加 而 逐渐 增加 ， 氢 化物 密度 为 2311 ugg 时 ， 模 型 届 服 应 力 为 8.23 GPa， 届 服 应 变 为 0.086， 杨 氏 模 量 
为 115.08 GPa。 

采用 COMB3 势 函数 计算 的 氧化 物 密度 对 Zr 基体 应 力 -应 变 曲 线 的 影响 ,与 AL Lloydi3 和 Ghaffarian U^! 
采用 MEAM 势 函 数 计算 的 氧 浓度 对 Zr 力学 性 能 的 影响 趋势 一 致 ; 对 纯 Zr 模型 计算 得 到 的 届 服 应 力 (13.31 
GPa) 5 Amit KH EAM 势 函数 计算 得 到 的 数值 (14 GPa) 相 比 误差 在 5% 之 内 , 说明 COMB3 势 也 适合 
用 来 进行 错 合 金 力学 性 能 的 计算 。 


—141— Pure Zr 
—$— 154 ug/g 
—4— 308 ug/g 
—v— 462 ug/g 
—*— 616 ug/g 
ーー 770 ug/g 
—>— 924 ug/g 
—*— 1078 ug/g 
ーー メー 1232 ug/g 
—*$— 1386 ug/g 
—3$9— 1540 ug/g 
ーー 1695 ug/g 
ーー プー 1849 ug/g 
ーー 2003 ug/g 
ーー 2311 ng/g 


Stress/ GPa 


0.00 0.05 0.10 0.15 0.20 0.25 0.30 


Strain 


图 3 不同 気化 物 密度 的 Zr 基体 所 对 应 的 应 力 -应 变 曲线 


Fig.3 Stress-strain curves corresponding to Zr substrates with different hydride densities 
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一 = 一 Yield stress 
— ョ ーー Young’ s modulus 


Yield Stress/ GPa 
Young' s modulus/ GPa 


ウシ ンー ーー キーー ョ ーーー ョ ーー ドー - 


6 
Pure Zr 154 308 462 616 770 2924 1078 1232 1386 1540 1695 1849 2003 2157 2311 
Hydride density / ug:g' 


图 4 氢化 物 密度 与 屈服 强度 和 杨 氏 模 量 的 关系 
Fig.4 The yield strength and Young's modulus of the Zr-based model correspond to different hydride densities 


2.2 ”微观 结构 分 析 


图 5 显示 了 不 同 密度 的 氧化 物 沉淀 (154~770 ugg) 拉 伸 模拟 中 的 微观 结构 演变 及 应 力 云图 。 由 图 中 
微观 结构 演变 可 以 看 出 ， 位 错 缺陷 均 沿 着 氧化 物 沉 演 界 面 形 核 ， 并 且 随 着 氧化 物 密度 的 增加 ， 氢 化 物 与 Zr 
基体 的 接触 面积 也 在 增加 ， 随 之 产生 的 缺陷 也 越 多 。 随 着 氧化 物 密度 的 增加 ， 位 错 越 来 越 多 ， 并 且 位 错 大 
多 是 单位 错 ， 使 得 模型 拉 伸 强 度 不 断 下 降 。 在 氧化 物 密 度 为 1078 ng/g 时 出 现 转折 ， 位 错 的 大 量 增 加 ， 使 
得 其 缠绕 在 一 起 ， 出 现 了 位 错 塞 积 现象 ， 模 型 的 拉 伸 强度 开始 增加 。 由 于 塞 积 并 不 剧烈 ， 所 以 模型 力学 参 
数 只 是 小 幅 增加 。 由 图 中 的 应 力 云 图 可 以 看 出 ， 随 着 氧化 物 密度 的 增加 ， 氢 化 物 周围 产生 的 应 力 集中 区 域 
也 在 扩大 ， 这 也 是 由 于 在 氧化 物 周 围 产生 位 错 导致 出 现 应 力 集中 ， 位 错 密度 的 不 断 增 加 使 得 扩大 应 力 集 中 
区 域 


ur 
is 


为 了 研究 不 同 密度 氧化 物 在 变形 机 制 中 的 作用 ， 对 体系 } 行 拉 伸 时 氢化 物 端 部 应 力 值 和 位 错 线 长 度 变 
化 统计 如 图 6 和 图 7 所 示 。 由 图 6 可 知 ， 随 着 氧化 物 密度 的 增加 ， 同 一 应 变 值 下 氧化 物 端 部 应 力 值 随 之 增 
加 ， 氧 化 物 密度 为 154 ug/g 时 ， 気化 物 3 省 部 应 力 信 最 大 为 5.899 GPa, 氧化 物 密度 为 770 pgg 时 ， 氧 化 物 
We o 这 与 图 5 中 気化 物 密度 増加 , 氨 化 物 周围 应 力 集中 加 剧 , 应 力 集中 区 域 扩 
大 的 微观 演变 一 致 。 由 图 7 可 知 ， 模 型 出 现 位 错 所 对 应 的 应 变 值 随 着 氧化 物 密度 的 增加 逐渐 变 小 。 毛 化物 
$0 UR 位 错 在 应 变 为 15.5% 时 开始 产生 并 增长 ， 相 较 于 届 服 点 (应 变 值 为 14.75% ) 
提前 5.1%， 位 错 线 密度 峰值 为 1478.9287 Å; 氧化 物 密度 为 154 ugg 时 ， 位 错 在 应 变 为 11% 时 开始 产生 并 
增长 ， 相 较 于 屈服 点 〈 应 变 值 为 12.55%) 提前 12.35%， 位 错 线 密度 峰值 为 5532.16308 À; 氢化 物 密度 为 
770 ugg 时 ， 位 错 在 应 变 为 4.5% 时 开始 产生 并 增长 ， 相 较 于 届 服 点 (应 变 值 为 9.4%) 提前 52.12%， 位 错 
线 密度 峰值 为 8377.27 A。 这 表明 随 着 氢化 物 密度 的 增加 , 位 错 出 现 的 应 变 点 将 随 之 提前 ， 并 且 位 错 线 密度 
随 氧化 物 密度 的 增加 而 升 高 ， 这 与 图 5 的 微观 结构 演变 一 致 。 
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图 5 不 同 氢化 物 密度 所 对 应 拉 伸 模拟 过 程 中 的 微观 结构 演变 , 氧化 物 密度 为 (a) 154 ug/g» (b) 308 ug/g, (c) 402 ug/g, (d) 
616 ug/g, (e) 770 ug/g 
Fig.S Microstructure evolution during tensile simulation with different hydride densities 


(a) 154 ug/g, (b) 308 ug/g, (c) 402 ug/g, (d) 616 ug/g, and (e) 770 ug/g 
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Hydride end stress value/ GPa 


Fig.6 
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图 6 不 同 氨 化 物 密度 的 拉 伸 模拟 中 氧化 物 尖 端 应 力 值 


Hydride tp stress values in tensile simulation with different hydride densities 
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图 7 不 同 氢化 物 密度 的 拉 伸 模拟 中 位 错 线 长 度 变化 曲线 


Fig.7 Length variation curve of dislocation line in tensile simulation with different hydride densities 


3 ”结论 


本 工作 通过 构建 含 氧化 物 的 模型 ， 采 用 COMB3 势 函数 ， 利 月 


模拟 ,探究 了 氧化 物 密度 对 力学 性 能 的 影响 ， 得 到 以 下 结论 : 
(1) 随 着 氢化 物 密度 的 增加 ,Zr 基体 的 屈服 强度 、 应 变 和 杨 氏 模 量 降低 ， 在 氧化 物 密度 增加 至 1078 


ug/g b], 各 力 学 情 


分 子 动力 学 方法 对 Zr 模型 进行 单 轴 拉 伸 


FE 能 指标 降 至 最 低 ， 之 后 ， 随 氧化 物 密度 增加 ， 各 力学 性 能 指标 逐渐 升 高 ，; 


(2) 気化 物 密度 的 増加 , 使 気 化物 与 Zr 基体 接触 面积 增 大 ,使 应 力 集中 增 大 ， 有 利于 位 错 形 核 ， 使 


Zr 基体 的 屈服 强度 与 届 服 应 变 下 降 ， 当 位 错过 多 时 ， 位 错 缠绕 造成 塞 积 ， 导 致 基体 的 强度 升 高 。 
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